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We present recent results from the BABAR collaboration in several areas of research. These 
include searches for new physics via measurements of radiative-penguin B decays, lepton- 
number and lepton-flavor violations in B decays, and CP violation in tau lepton decays. 

1 Introduction 

Until 2008, the BABAR experiment recorded e + e~ collisions at the T(4S) resonance with 471 x 10 6 
BB pairs produced (corresponding to an integrated luminosity of 429 fb _1 ). Four years after 
the end of the data-taking period, BABAR is still producing many new and interesting results. 
Despite the start of the LHCb experiment at CERN, BABAR is still competitive, especially for 
channels involving neutral particles (such as photons, 7r° or Kg), tau leptons, and neutrinos. 

We present here a selection of recent results from the BABAR experiment. We show exclusive 
measurements of the b — > sj transition rate as well as a study of the photon energy spectrum. 
Then we present an analysis of the angular distributions in the decay B — > K*l + l~. We detail 
a search for lepton-number violating processes in B + — > h~£ + l + decays, as well as a search for 
the decay B^ — > h^ri. Finally, a search for CP violation in the t~ — > ir~ ~Kg (> 07T°) v T channel 
is presented. 

2 Exclusive measurements of b — > s'j transition rate and photon energy spectrum 

As it is well known, flavor changing neutral currents (FCNC), such as b — > sj, are forbidden 
at tree level in the standard model (SM). However, FCNC of this type are predicted to occur 
at loop level with the following rate^4 B(B — > X s j) = (3.15 ± 0.23) x 10 -4 , for a minimum 
photon energy E~ > 1.6 GeV measured in the B meson rest frame, and where X s is the final 
state of the s quark hadronic system. The world average experimental valued is measured at 
B(B X s j) = (3.55 ±0.25 ±0.09) x 10~ 4 , for £ 7 > 1.6 GeV, and where the second uncertainty 
is due to the extrapolation from the experimental photon energy (between 1.7 and 2.0 GeV 
depending on the experiments) to 1.6 GeV. The calculation is performed at next-to-next-leading 
order in the perturbative term, with the first order radiative penguin diagram for the b — > S7 
transition having a W boson and a t, c, or u quark in the loop. Comparing the experimental 
and predicted values allows a precision test of the SM. Particles from new physics could enter in 
the loop, and would affect this transition rate. Furthermore, the photon energy spectrum from 
this reaction gives insight into the momentum distribution function of the b quark inside the 



B meson, and helps to constrain the uncertainty on the Cabibbo-Kobayashi-Maskawa matrix 
element V u b- 

In this BABAR analysis^, we use a "sum of exclusives" approach, where we reconstruct 38 
different X s final states (listed in Table [1]) . By fully reconstructing X s , we obtain the energy of 

the transition photon in the B rest frame with = \ m ^ a ■> where is the B mass, and 
mx s is the invariant mass of the X s system. We use a range of 0.6 < mx a < 2.8 GeV/c 2 , which 
corresponds to a photon energy range of 1.9 < < 2.61 GeV. In the analysis, this range is 
divided in 18 different mx a regions. 

Two types of signal Monte Carlo (MC) events are generated: one for the K*(892) region 
(corresponding to mx a < 1-1 GeV/c 2 ) and one for the inclusive region (corresponding to 1.1 < 
mx a < 2.8 GeV/ c 2 ). A flat photon spectrum level is used in the inclusive region at the generation 
level, which allows us to reweight to match whichever spectrum model we choose. 

We use three classifiers (based on random forest classifiers) to reject the background, each 
optimized in four different mx s regions. These three classifiers help respectively to choose the 
best B candidate, to veto photons coming from a neutral pion, and to reject the continuum 
background. The signal yield is extracted in each mx s region from a fit to the beam-energy 

substituted mass, uies = \J (\/^/2) 2 — {P*b) 2 i where p* B is the momentum of the reconstructed 
B meson in the center-of-mass. An example of such a fit is shown in Fig. [U 

The fragmentation of the hadronic system is modeled with JETSET with a phase-space 
hadronization model. We observe some differences between fragmentation in the MC and data 
samples. To correct this, we group the final states by topology and correct the signal contribu- 
tion in the MC to better reflect the data. Furthermore, several quark hadronization models are 
tested and the observed differences are included in the systematic uncertainties. 

The result of the measurements of the partial branching fractions in each mx a region is 
shown in Fig. CD We fit this spectrum with two models: the kinetic model and the shape- 
function model : The result is shown in Fig. [1] for the kinetic model. We have also extracted 
the parameters from the heavy quark effective theory (HQET) for these two models, obtaining 
values compatible with the world average. The mean and variance of the photon energy spectrum 
have also been extracted, which can be used to constrain different kinds of models. Summing 
the 18 bins, we find the total branching fraction for E^ > 1.9 GeV to be B{B —> X s j) = 
(3.29±0.19±0.48) x 10 -4 , where the first uncertainty is statistical and the second is systematic. 



Table 1: The 38 final states of the B meson reconstructed in this analysis. 
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3 Angular distribution in the decays B — > K*i + l 

The decay B — > K*£ + £~ occurs in the SM via penguin (electromagnetic penguin loop or elec- 
troweak Z° penguin loop) and W + W~ box diagrams. The effective Hamiltonian of these re- 
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Figure 1: Left: example of tues fit for the region 1.4 < mx s < 1.5GeV/c 2 . Points with statistical errors are the 
data, the solid line represents the total fit, the thick-dashed line represents the signal, and the other lines show 
background contributions. Right: partial branching fractions in each mx s region (dots) and fit of the spectrum 

using the kinetic model (solid line). 



actions depends in particular on the Wilson coefficients Cj. Particles from new physics could 
enter in the loop (such as charged Higgs, squarks, neutralinos, and charginos) and could lead to 
sizeable deviations from the SM predictions. In this section, we focus on an analysis^ on the 
angular observables in B — > K*£ + £~ (the branching fraction and rate asymmetries are addressed 
elsewhere in these proceedings^). The angular distribution depends in particular on the angle 
6k, the angle between the K and the B in the K* rest frame, and on 9g, the angle between 
the £ + (£~) and the B (B) in the £ + £~ frame. From these angles, we can obtain the fraction of 
the longitudinal polarization of the K*, Fl, and the lepton forward-backward asymmetry, Afb- 
The predicted distributions depending on these quantities, as a function of the dilepton mass 
squared s = rrij +e _ , read: 

iy^p = lF L coJe K + -Al-F L ){l-co S 2 e K ), 
r a cos Ok 2 4 

^T^- = ~ cos2 Ot) + l(l-F L )(l + cos 2 e ) + A FB cos 9 e . 

1 a cos 0£ 4 8 

In the SM, at low s, where the effective Wilson coefficient Cf dominates, Afb is expected to 
be small, crossing the zero axis around s ~ 4 GeV 2 /c 4 . At high s, the product of C| ff and Cfg 
is expected to give a large positive asymmetry. Contributions from new physics could change 
dramatically this picture. 

We use seven bins in the s variable, with the same binning as the experiments CDF, Belle, 
and LHCb, in order to ease the comparison and the average. Two regions in s are dominated 
by J/ip and tp(2S) and are used as control samples (and are vetoed in the analysis). We use 
the decays B + -> K*+l + l~, followed by K* + -> K + ir , K®-k + , and B° -> K*°£ + £~, followed 
by K*° -»• K+-K-, with £ = e,/J,. The decay B + K* + fi + fi-, K* + -)• K + tt° is not used since 
it showed no improvement to the analysis. The continuum and BB backgrounds are rejected 
using bagged decision trees (BDT), based on AE (the difference between the expected B energy 
and the reconstructed B energy), event shape and vertexing variables. A likelihood ratio R is 
constructed from the BB BDT. The angular observables are extracted from simultaneous fits 
over the combinations of final states. The strategy in each bin in s is three-fold: first, for each 
of the five decay modes, we fit the variables thes, M(Ktt) (the reconstructed invariant mass of 
the K*), and R, where we fix the fit parameters for the next stage; second, we fit the cos 9k 
distribution to extract Fl) and third, we fix Fl, and fit the cos#£ distribution to extract Afb- 

The results are presented in Fig. We obtain significantly more precise values than the 
ones from either Belled or CDfE. Other experiments, especially LHCb^, are dominated by the 



B° — > K*°£ + £ channels. On the contrary, we are able to study B + — > K* + £ + £ and notice 
some tension at low s in this channel. 




Figure 2: Distributions of Fl (left) and Afb (right) in the seven bins in the s = rnf +e _ variable. The two 
vetoed regions correspond to regions dominated by the J/ijj and tp(2S) resonances. The blue downward triangles 
correspond to B + — > K* + £ + £~ , the red upward triangles to B° — > K*°£ + £~ , and the black dots to the combination 
B — > K*£ + £~ . The blue dotted line shows the SM prediction (the theoretical uncertainties are of the order of 

5-10% at low s and 10-15% at high s). 



4 Search for lepton-number violating processes in B + — > h £ + £ + decays 

In the SM, the lepton number L is conserved in low-energy collisions and decays. However, 
since the neutrinos are oscillating, we know that these particles have a mass. If the neutrinos 
are of Majorana type, then L violation becomes possible. This violation could be seen in B 
meson decays via the process B + — > h^£ + £ + (see Fig. [3]). This process becomes resonant if 
the neutrino mass lies between the h meson and the B meson masses. Such processes involving 
meson decays ^ are an alternative to neutrinoless double beta decays. Beyond the SM, this 
diagram could lead to a lepton violation via new physics, such as left-right symmetric gauge 
theories, SO(10) supersymmetry, i?-parity violating models, or extra-dimensions. 




Figure 3: An example of a diagram of the process B + — > h £ + £ + with lepton-number violation. 

In this analysis four final states are considered: the decay channel is B + — > h~£ + £ + 
with h = K, 7r and I = e, \x. The selection is identical to the previous section since the final 
states are very similar. Unbinned maximum likelihood fits of twes an d R, the likelihood ratio, 
are performed for each of the four modes. We use the B + — > J/iph + data control sample to 
obtain the mES fit parameters. The fits are shown in Fig. [4j where we observe that no signal 
is seen. From these negative results, we obtain upper limits at 90% confidence level (CL) 
on the four channels: B{B + -> K~e + e + ) < 3.0 x 10~ 8 , B(B + -> R-^ + n + ) < 6.7 x 10~ 8 , 
B(B + -»■ 7r-e+e+) < 2.3 x 10~ 8 , and B(B+ -»■ TT-fi + fi + ) < 10.7 x 10" 8 . These upper limits are 
40-70 times more stringent than previous limits set by other experiments. These results can be 
translated on upper limits on the branching fractions as a function of the mass rri£+f l - (Fig. H]), 
which can be related to the Majorana neutrino mass for diagrams of the type of Fig. [3l 




Figure 4: Left: jties distributions for a) B + — > K~e + e + , b) B + —> K~/j, + fi + , c) B + — > 7r~e + e + , and d) 
B + — > tt~ (i + )i + . The solid blue line is the total fit, the dotted magenta line is the background, the solid green 
histogram is the signal. Right: 90% CL upper limits on the branching fractions as a function of the mass m e + h - . 



5 A search for the decay modes — > h iz r£ 

FCNC and charged lepton flavor violation are forbidden in the SM at tree level. However, in 
many extensions of the SM, these effects could be enhanced, especially for the second and third 
generation We study the eight final states — > h ziz r£, with h = K,tt and £ = e,[i. 
The final states B^ — > K^re, B^ — > tt^t^l, and — > 7r ± re have never been done before. 
We fully reconstruct the hadronic B on one side (the "tag" B) using final states of the type 
B~ —7- D^°X~ , where X~ is composed of 7r ± , K^, K®, and it . This determines the three- 
momentum of the other B (the "signal" B) on the other side and thus allows us to indirectly 
reconstruct the r lepton through: 

Pt = -Ptag - Ph - Pi, 
E T = -Ebeam ~ Eh — E(_, 

m T = \/E%- \p T \ 2 , 

where (E T , p T ), (Eh, Ph), and (En, pi) are the corresponding four-momenta of the reconstructed 
signal objects, and where ptag is the three-momentum of the tag B, and E^ e£im the beam energy. 
The r is required to decay to a "one-prong" final state: r — > evV, r — > iivV, and r — > (mr )^^ 
with n > 0. The signal branching fraction is determined by using the ratio of the number of 
B^ — > h^ri signal candidates to the yield of control samples of B + — > D^°£ + i'; D° — > K + tt~ 
events from a fully reconstructed hadronic B^ decay sample. 

The background is mainly coming from semileptonic B decays when the charge of the signal 
B is the same as the one from the primary lepton, and mainly from semileptonic D decays when 
the charges are opposite. We remove these backgrounds by rejecting the signal B candidates 
where two of their daughters are kinematically compatible with originating from a charm decay. 
After this requirement, we reject the continuum background using a cut on the likelihood ratio 
R, based on particle identification and event shape variables. 

The signal region is defined as ±60MeV/c 2 around the indirectly reconstructed r mass m T . 
Figures [5] and [6] show the m T distributions for the data, for the background, and for the signal 
MC. No signal is observed, which allows us to put 90% CL limit on the branching fractions. 
Assuming B(B + — > h + r~£ + ) = B(B + — > h + r + £~), we obtain the combined limits shown in 
Table [2j These limits can be translated into model-independent bounds on the energy scale of 
new physics in flavor-changing operators ^: Ag d > 11 TeV and A^ s > 15 TeV (at 90% CL), 
which improved the previous limits of 2.2 and 2.6 TeV, respectively. 
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Figure 5: Observed distributions of the r invariant mass for the B —¥ KtI modes. The distributions show the 
sum of the three r channels (e, 7r). The points with error bars are the data. The solid line is the background 
MC which has been normalized to the area of the data distribution. The dashed vertical lines indicate the m T 
signal window range. The inset shows the m T distribution for signal MC. 




Figure 6: Observed distributions of the r invariant mass for the B — > nr£ modes. The distributions show the 
sum of the three r channels (e, /i, n). The points with error bars are the data. The solid line is the background 
MC which has been normalized to the area of the data distribution. The dashed vertical lines indicate the m T 
signal window range. The inset shows the m T distribution for signal MC. 



Table 2: Branching fraction central values and 90% CL upper limits (UL) for the combination B(B + — > h + r£) = 



-j. h + T~£ + ) + g(B + ->■ h + T + £~). 
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6 Search for CP violation in the decays r — > ir (> 07r°) v T 

CP violation, until now, has only been observed in hadronic decays {K, B, and D systems). 
However, Bigi and Sanda predict^] a non-zero decay rate asymmetry for r decays to final states 
containing a Kg meson, due to the CP violation in the kaon sector. The decay rate asymmetry 
is: 

_ r (r+ -> it + K° V T ) - r (t - -> vr-K° zy r ) 
74(3 = r (r+ vr+K° F T ) + r (r- ->• vr-K° i/ T ) ' 

and is predicted to be equal to (0.33 ± 0.01)%. Any deviation from the SM prediction could be 
a sign of new physics. It has to be noted that A q is independent of the number of neutral pions 
in the final state. Recently, Grossman and Nir^^ noticed that the calculation needs to take 
into account interferences between the amplitudes of intermediate K® and K\ mesons, which 
are as important as the pure K® amplitude. This means that Aq depends on the reconstruction 
efficiency as a function of the K® — > 7r + 7r~ decay time. 

We study here^^the decay channel r~ — > ir~Kg (> 07r°) u t . The event is divided into two 
hemispheres, one corresponding to the signal side, and one to the tag side with t~ — > £~U£V T , 
£ = e,fj,. The selection of the signal events requires that the invariant mass of the reconstructed 
r lepton is smaller than 1.8GeV/c 2 . A first likelihood ratio is constructed to reject continuum 
background based on energy, calorimeter clusters, thrust, and momentum. A second likelihood 
ratio is aimed at reducing the K® background based on K® reconstruction parameters. After 
the selection, we obtain 199064 candidates for the electron tag channel (e-tag), and 140602 for 
the muon tag channel (ii-tag). The composition of the sample in term of signal and background 
events is presented in Table El 



Table 3: Composition of the sample after all selection criteria have been applied. 



Source 


Fractions (%) 
e-tag A*-tag 


T~ -^tt-^(> 0vr u )i/ T 
T~ -»■ K~ K° s {> Ovr )^ 
T~ -> TP" K°K° u T 
Other background 


78.7 ±4.0 
4.2 ±0.3 
15.7 ±3.7 
1.40 ±0.06 


78.4 ±4.0 
4.1 ±0.3 
15.9 ±3.7 
1.55 ±0.07 



We need to correct the raw asymmetry from the pollution of the other modes shown in 
Table El For the mode t~ — > K~ K®{> Qit®)v T , the expected asymmetry is opposite to the 
one from the signal, and for the mode r~ — > ir~ K°K° u T , the expected asymmetry is zero. 
Furthermore, an additional correction was pointed out recently^l we need to take into account 
a correction on the asymmetry Aq due to the different nuclear-interaction cross-section of the 
K° and mesons with the material in the detector. We calculate this correction to be (0.07 ± 



0.01)% and we subtract it from the measured asymmetry. After all corrections are applied, and 
after combining the results from the e-tag and //-tag, we obtain Aq = (—0.36 ± 0.23 ± 0.11)%. 

As we mentioned, this result should be compared with the prediction of the SM, corrected 
by the K® — > tt + tt~ decay time dependence. Using the MC sample, we find a multiplicative 
factor of 1.08 ± 0.01. The SM decay-rate asymmetry, after correction, is then predicted to be 
A| M = (0.36 ± 0.01)%. We observe that our measurement is 2.8 standard deviations away from 
the SM prediction. 

7 Conclusions 

The BABAR collaboration, four years after its shutdown, is still producing competitive results, 
with many more to come in the near future. We have presented here five new results. In general, 
we have a good agreement with the prediction of the SM, except for two studies: we observe 
some tensions at low wn%+t- for Fl and Afb in the channel B + — > K* + £ + l~; and we measure 
a CP violation parameter 2.8 standard deviations away from the SM prediction in the channel 
r~ — > n~K® (> 07T°) v T . The actual statistics is not sufficient to tell whether or not these could 
be indication for new physics. 
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